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Lattice and continuum fluid models with repulsive-core interactions typically display a domi- 
nant, critical-type singularity on the real, negative activity axis. Lai and Fisher recently suggested, 
mainly on numerical grounds, that this repulsive-core singularity is universal and in the same class as 
the Yang-Lee edge singularities, which arise above criticality at complex activities with positive real 
part. A general analytic demonstration of this identification is presented here using a field-theory 
approach with separate representation of the repulsive and attractive parts of the pair interactions. 
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I. INTRODUCTION 

In 1952 Lee and Yang Q| proved that the zeroes of 
the partition function of a ferromagnetic Ising model (or 
equivalently, a lattice gas with single-site hard cores) as a 
function of complex magnetic field H (or chemical poten- 
tial /i) are confined to the imaginary H axis for real tem- 
peratures T. Yang and Lee argued that for a system 
above its critical temperature Tc, the partition function 
must be nonzero throughout some neighborhood of the 
real axis in the complex reduced magnetic-field plane. 



H/ksT = h' + ih". 



(1) 



Thus for T > Tc, a gap free of zeros will be found on the 
imaginary h axis with edges at, say, ±ih^{T). Equiva- 
lently, in fluid-language, for T > Tc, there will be a gap 
in the complex activity plane with edges at 



(2) 



as illustrated in Fig. 1. Here we define the reduced ac- 
tivity in d spatial dimensions by 



Wo /A 



d 



(3) 



where (3 = l/kgT, while vq is a microscopic reference 
volume (taken as the cell volume for a lattice gas) and 
is the thermal de Broglie wavelength. 
If one defines the density of Yang-Lee zeros, g{h"), 
so that when N, the number of spins (or lattice sites) 
becomes infinite, N g{h")dh" approaches the number 
of zeros between ih" and i{h" -I- dh") on the imagi- 
nary h axis, one must have g{h") = for \h"\ < 
K{T), for T > Tc . Kortman and Griffiths |] 
pointed out that the density of zeros beyond the gap 
should be expected to exhibit a power law singularity — 
the Yang-Lee edge singularity — of the form 

g{h")r.\\h"\-hATr for \h"\^K{T)-. (4) 

This singularity proves analogous to an ordinary criti- 
cal point obeying scaling laws and exponent relations 
although there is only one relevant scaling field [3, 4]. 



Thus the basic exponent cr in (4) satisfies a = 1/6 = 
{d~2 + r^)/{d + 2--rj) §. 

In a general renormalization group analysis, it was 
shown 0] that the field theory controlling the Yang-Lee 
fixed point is described by a pure imaginary iip^ cou- 
pling. This leads to a critical dimension dc — 6 above 
which the classical, mean-field value c = ^ applies. 
To first order ine — 6 — d>0 one has a{d) 
and T] 
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ge. Expansions of a{d) to order e are known 
and one has a{l) — — ^ Jp|, o-(2) = — |, and finds, 
numerically, a{3) ~ 0.088 : see Lai and Fisher @] 
who review previously known relationships of the Yang- 
Lee edge singularity to a number of different problems, 
specifically, isotropic branched polymers or, equivalently, 
undirected lattice animals with or without loops allowed, 
in (d + 2) dimensions [8, 9]; Anderson localization [ p^ ; 
and directed branched polymers (or directed, loop-free 
lattice animals) in d -I- 1 dimensions pl| . 

By contrast, an apparently quite different type of sin- 
gularity arises in fluid systems when the particle inter- 
actions have repulsive cores. If the pair interaction po- 
tentials are purely repulsive (i.e., positive), the (reduced) 
cluster integrals, 6„(T), in the activity or fugacity series 
for the (reduced) pressure, namely. 



- 2^bn{T)z 

n=l 



(5) 



are known ||T^ to alternate in sign: this implies a domi- 
nant singularity on the negative z-axis that determines 
the radius of convergence, R, of the series: see Fig. 1. In 
the vicinity of this singularity, say at z = zq{T) = —R, 
the reduced pressure can be written as 



p{z) = Pq + pi{z - Zq) + P2{z ~ Zq^ + ■ ■ 

+ Piz- zo)^[l + agiz- zo)" + 



(6) 



where the exponents (p and 9 are anticipated to be non- 
integral. 

Indeed, in 1984, Poland [|l3| studied a variety of lat- 
tice models and a continuum fluid of hard squares and 



proposed that this repulsive-core singularity is charac- 
terized by a universal exponent <^(rf). Subsequent con- 
firmation came from Baram and Luban who inves- 
tigated further models including dimers on lattices, par- 
allel hypercubes in continuum space and the soft-core 
single-component Gaussian-molecule model. More re- 
cently, Lai and Fisher |Q found similar behavior in a 
binary Gaussian-molecule mixture using very long series 
expansions for d = 1, ■ • ■ , 6. In that case the singularity 
at z = zo was drawn out into a continuous locus; but the 
estimates for (j){d) supported the universality hypothe- 
sis. (Precise estimates, for the leading "correction-to- 
scaling" exponent 6{d), in (6) for all d, including ^1) = 1, 
0{2) = |, and e{3) ~ 0.62, were also generated 0.) 

However, Lai and Fisher noticed in particular 
that, when compared with previous knowledge about the 
Yang-Lee edge exponent g(c?), the exact results for 
for d = 1 and d = oo and for d = 2 and 
the various numerical estimates j?) for d > 2, strongly 
suggested the identification 

(j){d) = a{d) + 1 for all d. (7) 

Hence they proposed that the universal repulsive-core 
singularity belongs to the Yang-Lee edge critical univer- 
sality class. 

In analytical support of this identification (for general 
d) they appealed to earlier work by Kurtze and Fisher 
[ pTf who had proved that when T ^ oo the Yang-Lee 
edge singularity in ferromagnetic Ising models precisely 
describes the dominant singularity on the negative z axis 
of a fluid of hard dimers on the same lattice, each dimer 
occupying one bond and two adjacent lattice sites. This 
correspondence with a dimer gas was carried further by 
Shapir who used field-theoretic arguments to demon- 
strate the identity generically for all T > Tc- However, 
dimers are rather special objects with orientational de- 
grees of freedom (although these seem, post facto, to 
have no effect on (j>{d) jlj]). Furthermore, Shapir's field- 
theoretic approach was rather special and did not seem 
extendable to more general particles with soft repulsive 
cores, with the additional presence of attractive forces, or 
to systems also displaying critical behavior for positive z 
(as illustrated in Fig. 1). 

In this article we repair this gap in the theory. Specif- 
ically, we consider a general single-component fluid with 
a pair interaction potential, U{r), which contains both 
repulsive and attractive parts: of course, repulsive terms 
are always essential to ensure thermodynamic stability. 
To be concrete and explicit we analyze lattice systems 
in which multiple occupancy of a site is forbidden; how- 
ever, it must be stressed that the repulsive interactions 
we consider are not confined to such trivial single-site 
hard cores. On the contrary, essentially we suppose only 
that the potential U{v) is of finite range and may be de- 
composed according to 

;7(k) = (^'""^UiR) = i&(k) - vik), (8) 

R#0 



where the sum runs over lattice sites R = r-,- while the re- 
pulsive and attractive parts, w(k) and w(k), respectively, 
are both positive (if they do not vanish identically). More 
specifically we will use 

w{k)=wo{l-k^a^ + ---)>0, (9) 



v{k) ^ vo{l - + ■■■)> 0, (10) 

so that a and b represent the interaction ranges of the 
repulsive and attractive components. The real-space po- 
tentials acting between sites j and k (j ^ k) , namely, 
Wjk = w{vk — Vj) and Vjk = v(yk — fj), follow by Fourier 
inversion. For d = 3 one may imagine Yukawa forms: 
w{v) = Woe~''/''/r and v{t) = Vbe~''/''/r' with a <b\ but 
that is certainly not essential. Likewise, the leading lat- 
tice isotropy assumed for convenience in (9) and (10) is 
not necessary. Our treatment extends straightforwardly 
to multicomponent fluids and, at least formally, gen- 
eralizes readily to continuum systems. 

On the basis of (8)-(10), we develop a field-theoretic 
analysis and show that there is, in general, a repulsive- 
core singularity at some zq{T) on the negative activity 
axis (see Fig. 1) that is described within an LGW renor- 
malization group framework by a fixed point Hamiltonian 
with a purely imaginary cubic coupling, and, hence, 
lies in the same universality class as Yang-Lee edge crit- 
icality. A more or less novel aspect of our treatment 
is that separate — sinc-Gordan [|9| and Kac-Hubbard- 
Stratonovich (KHS) ||2C|-[2^ — transformations are used 
for handling the repulsive and attractive parts of the in- 
teraction potential: compare with Refs. 23 and 24. 



II. FORMULATION 



For a lattice of volume V with sites labeled j,k,l — 
1, 2, • • •, TV = V/vo, let pj = or 1 according as site j 
is or is not occupied by a particle. Then, recalling (3) 
and (8)-(10) et seq., the grand partition function for the 
lattice system is 



Z{T,z) = Tr^ <^ z^^'^ exp[-i/?^p,(«;,fe - Vjk)pk] \ ■ 

(11) 

Now, utilizing the positivity of w{k) and v{k), we may 
apply a sine-Gordon transformation |19| to the repulsive 
terms, Wjk, and a KHS transformation ]2C|-[24| to the at- 
tractive terms Vjk (where, as usual, it is most convenient 
to utilize periodic lattice boundary conditions) . Neglect- 
ing an unimportant constant factor, this yields 



2 



z = 



f Vx 



l.m 



where w ~ [[3wjk\ — [wj/}]^^ and similarly for v. Per- 
forming the trace over the pj then yields the transformed 
reduced Hamiltonian 

T-^iv^ X] = 5 J2^'f'i^Jk'Pk + XjVjkXk) (13) 
-^ln(l + ze-*'^^+'^0- 



Now, in seeking to understand the possible singular- 
ities in the reduced pressure p{T,z) = {vo/V)lnZ, we 
will, initially, neglect fluctuations and study the saddle 
point(s) which extremize the integrand exp{— 7i[(/?, x]} 
in (12). We expect spatially uniform solutions ipj — ipo, 
Xj — Xo (all j) to sufhce: from dTL/dipj = one thus 
finds 



ze 



1 + ze^*'^"+xo 



0, 



(14) 



for the repulsive terms. On premultiplying by the ma- 
trix w = [I3wim] (see, e.g., [25]) and using the Fourier 
representation (9), one obtains the simpler form 



Similarly from dJi/dxj = we find 

XO =/3«o/(l + 2-'e-(>^«-^'^«)). 



(15) 



(16) 



for the attractive terms. 

Before analyzing these coupled saddle-point equations, 
note that if one puts 



(17) 



the two equations combine simply and can be rewritten 
as 



^e-^ = -z{ij + l3Uo), 



(18) 



where the total strength of the pair potential is measured 
by 



Uo = U{0) = i&o ~ va- 



(19) 



Furthermore, by substitution of any solution of (18) on 
the right hand sides of (15) and (16) one obtains the sep- 
arate solutions ipo, which may evidently be imaginary, 
and Xo- 



iIII. TRIVIAL HARD CORES AND GAS-LIQUID 
CRITICALITY 



(12) We remark, first, that following the pioneering study 
of Hubbard and Schofield , all authors interested in 
ordinary gas-liquid criticality have treated the repulsive 
interactions in a fluid by use of a reference system: for 
a recent example, see Brilliantov Unless this refer- 
ence fluid is essentially trivial, as for single-site hard cores 
on a lattice, this entails increasingly detailed knowledge 
of the correlation functions of the repulsive-core system 
1^,^. Nevertheless, it will be helpful for us to make 
contact with this approach by, initially, neglecting the 
nontrivial repulsive interactions embodied in u;(r). If the 
attractions are also neglected the grand partition func- 
tion (11) reduces simply to Z = (1 + z)^ . This predicts 
a repulsive-core singularity at zo = ^1 (a-U T) with an 
exponent (j){d = 0) = O(log) Evidently, the actual 
dimensionality, d, plays no role. 

Now, following traditional treatments (e.g., [^) let 
us introduce attractive terms with ■O(k) > 0. The KHS 
transformation then leads to the saddle-point equation 
(18) in which, now, Uq = —Vq < is negative and we can 
identify -0 directly with xo (since wq = and J T>ip etc., 
can be ignored). Fig. 2 then provides a graphical rep- 
resentation of (18) which can be readily analyzed: note 
that the bold curve depicts tfje~'^ which has a point of 
inflection at ■0c = 2 and a corresponding tangent that 
intersects the axis at ip = 4: see the dashed line. 

By inspection, one then sees that when z increases from 
z = (along the real axis) for temperatures such that 
f3vQ < 4 there is always a single saddle-point solution, 
TpoiT, z), which varies analytically with T and z: see, e.g., 
the lines labelled (a) and (b) in Fig. 2. Conversely, for 
lower T, when (3vo > 4, there is a single analytic solution, 
TpoiT, z), for small z [as on the line (c)] but for an inter- 
mediate range of z three distinct solutions arise, as illus- 
trated by (d); finally, for larger z only a single solution 
remains: line (e). Evidently the three solutions merge at 
a bifurcation point determined by the inflection point at 
4>c = 2: this leads to PcVo = 4 and Zc = ~ 0.135. 
Expanding H[x] in powers of (/? — /?c), {z — Zc) and 
Sxj — Xj ~ Xc, all taken as real variables, shows that 
this saddle-point bifurcation simply represents the antic- 
ipated classical or mean-fleld gas-liquid critical point at 
/cbT° = jVq. As usual, the corresponding LGW Hamil- 
tonian can be used as a starting point in a fleld-theoretic 
renormalization-group (RG) treatment which then leads 
to all the standard results. 

On the other hand, for /3 < /3c, i.e., T > Tc, one can fol- 
low Ref. [Ql and discover two Yang-Lee edge singularities 
at complex z (with small imaginary parts when T is near 
Tc): at the saddle-point level these have cr = i; but the 
flxed-point Hamiltonian is controlled by an i{Sx)^ cou- 
pling leading, as explained above, to an RG e — 6 — d 
expansion [^. 

If, next, small repulsive terms, w > 0, are introduced. 
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the overall interaction parameter Uq = wq—vq, in (18) re- 
mains negative and the arguments proceed in essentially 
the same manner (although, in due course, the field (p for 
the repulsive terms would normally be integrated out). 
As expected, neither the usual gas-liquid nor the Yang- 
Lee edge singularities undergo any change in character. 
However, if wq becomes sufficiently large, Uq becomes 
positive and then, clearly, the previous analysis fails! In 
particular, as seen in Fig. 3 [lines (a) and (b)], for positive 
(real) 2, there is always only a single, smoothly varying 
saddle-point solution, ipo{T,z). Of course, this does not 
(necessarily) mean that the usual gas-liquid and Yang- 
Lee edge singularities are lost. Rather, the form of the 
Hamiltonian, 7i[(/5,x], represents an inadequate starting 
point for a perturbative RG approach: instead, it be- 
comes necessary to integrate out (at least to some degree) 
the repulsive terms, i.e., to perform some J Dip integrals 
and, thereby, make contact with the reference-fluid treat- 
ments |t25-E7f. 



IV. THE REPULSIVE-CORE SINGULARITY 

On the other hand, when the repulsions dominate, so 
that Uo > (as we will assume hereon), one sees from 
Fig. 3 [e.g., lines (c) and (d)] that for a real negative 
z there is a unique positive saddle-point solution, say 
^j-{T, z) < 1, that vanishes when 2 — > 0. However, when 
2 approaches zq{T) < (from above) a bifurcation point 
is reached at which the saddle-point must become com- 
plex: see the broken line in Fig. 3 that corresponds to 
2 = 2o(T). (At 2 = 2o a second, larger saddle-point solu- 
tion, '0+(2^j z), merges with ip-{T^ 2). For 2 much larger, 
additional real negative saddle-point solutions, tp'^, ap- 
pear but these are not relevant to the dominant repulsive- 
core singularity.) As evident from Fig. 3, the bifur- 
cation point is located by the tangent passing through 
tp = —(3Uo < 0: this leads to 



values (po and xo [following from (15)-(17)]. In terms of 
Sipj = (fij - (fio, 5xj = Xj - Xo, and Sipj = 5xj - iS(pj, the 
Hamiltonian truncated at fourth order becomes 



Ti-iSfySx] = 5 J^i'f^^jk'Pk + XjVjkXk) 



(22) 



where it is convenient to put 

0(T,2) = -2-ie-'^V(l + 2-ie-'^°)', 



(23) 



(which is real and positive for 2 ~ 20 when / = 0, 2, 4, • • •) 
so that ro — C2iT, 2) and 

50 = ^3(1 + Co), uo- 2lC4(l+4Co + Co)- (24) 
In the usual continuum approximation this becomes 

n = [ d'^xi^niSip^^) + ri25^(x)5x(x) + ir225x'(x) 



+ic^(VM' + ^r^xi^Sxf 
-5o<5V'^(x) +Mo<5V'^(x)], 



with = a? / (3wq and = b'^ / (3vq and 



(25) 



rii= ksT/wQ - C,2{z,T), 
r22= ksT/vo + C2{z,T), 
ri2= -K2{z,T). 



(26) 



Before proceeding generally, let us suppose that only 
repulsive terms act, i.e., vq = 0. Then we may drop 
the integrals / Vx and see, by (17), that Sip = —iSip; 
but note that ipoiT) remains real, positive and less than 
unity. The continuum Hamiltonian then reduces to 



< MT) = yUo{[l + (iksT/Uo)]'/^ - 1} < 1, (20) ^"[-^^J = J ^[2^o(M + 2C^(VM ^ ^9oSp + u^S^ ], 



and 



2o(r)--[l-^o(T)] exp[-MT)], 
«-l/e/3C/o[l + 0(l//3C/o)], 



(21) 



the last result applying for strong repulsions {pUo 
00). Clearly, this saddle-point bifurcation represents the 
repulsive-core singularity. 

Incidentally, if following Hauge and Hcmmer ||l5|] , we 
had treated the {d = l)-dimcnsional continuum hard- rod 
gas with additional infinite-range, infinitely weak repul- 
sive Kac potentials, we would, at this point, have found 
zq = — l/e/3[/o in precise accord with the exact (limiting) 
calculations |1^. But, of course, our saddle-point treat- 
ment is not restricted to d = 1 even though in the Kac 
limit it will also become exact. 

To complete the analysis we may now follow stan- 
dard procedures by expanding about the saddle point 



where the controlling coefficient is 

, , _ _ 1 A z\ [l~22oe^-/'o(^)] 
to{T, z) _ rn - ^ (^1 - - j ^ ^^^„(T)p ■ 



(27) 



(28) 



which is positive for small 2, decreases as 2 — ~\z\ in- 
creases in magnitude, and vanishes linearly when 2 — > 
zq{T)+ < 0. The fourth order coefficient, Uq, is positive 
and approaches 



Woo 



i(l-7Ao)(l-V'o + ^^o')/^o >0, 



(29) 



when 2 —>■ 2o(T), while the cubic term is purely imagi- 
nary with a coefficient approaching 



5oo = -g(l-V'o)(2-^o)/V'o <0. 



(30) 
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At this point we have thus reached the same stage as in 
the original treatment [Q of the Yang-Lee edge singulari- 
ties. The imaginary term i{6ip)^ dominates the behavior. 
A momentum-shell RG analysis reveals an upper critical 
dimension dc = 6, an exponent 77 « — for e = 6 — d > 0, 
and, using a unit cutoff in momentum space, a fixed point 
value of (7q « (e/SiK^y^'^ , where Kq is the area of the 
unit sphere at d = 6 dimensions. All other results follow 
as previously [4,7]. 

Returning now to the general case of (25), in which 
attractive interactions are present (i.e., vq > 0) but Uq 
remains positive, we anticipate that only one linear com- 
bination of 6ip and 5x will become critical. To check this 
at the saddle-point level, it suffices to put 



<5x(x) = (5x'(x) + ^(C2/r22)'5v?(x). 



(31) 



The quadratic (nongradient) terms in (25) then become 



(32) 



The coefficient r22{T, z) remains positive for all negative 
z but one finds 



t(r, z) = rii 



' 12 

?'22 



h{T,z) 



[1-f /3z)o(l-Vo)/V'o'] 



(33) 



for 2; ~ zq- thus by (28), t vanishes when z zq{T). 
Integration over the noncritical field (5x'(x) is hence jus- 
tified and may be carried out perturbatively. Apart from 
additive terms that are nonsingular near zq{T)^ this fi- 
nally leads to a renormalized LGW Hamiltonian 



—ih^Sif - igBiSipy + ■ ■ 



(34) 



where the fourth and higher order terms have been 
dropped. For vq sufficiently small, the renormalized cou- 
plings tR and gn will differ little from t and g and, in 
particular, ha and gji will be real, so that the i{S(py 
term once again dominates. 

Stability of the saddle point requires that the renormal- 
ized repulsive range, a^, in (34) be real (and positive). 
Roughly, we expect aj^ oc {a'^WQ — b'^vo): the positiv- 
ity of this factor then represents a previously unstated 
restriction on the initial potentials. In practice, how- 
ever, if this (or the more precise) condition fails, it may 
be sufficient, as explained in the discussion of gas-liquid 
criticality, to perform further, nonperturbative renormal- 
izations to dampen the attractive interactions for z < 
and so expose again a saddle-point representation of the 
repulsive-core singularity at zo{T). Certainly, we must 
expect the form (34) to represent the singularity when- 
ever it is actually realized in a system; and, then, clearly 
it must belong to the Yang-Lee edge universality class. 



V. CONCLUSIONS 

In summary, by using separate field-theoretic transfor- 
mations for the repulsive and attractive parts of the pair 
interactions in a fluid, we have demonstrated in general 
the presence of a universal repulsive-core singularity on 
the negative axis at a value zo{T) < 0: see Fig. I. The 
behavior of the pressure in the vicinity of zq is stated 
in (6). The singularity belongs to the same universal- 
ity class as Yang-Lee edge criticality, as proposed by Lai 
and Fisher [Q: see Eqs. (l)-(4) and Fig. 1. The basic 
critical exponents are related via (7) (see also 0); the 
borderline dimensionality is dc = 6; and an iip^ coupling 
characterizes the LGW fixed-point Hamiltonian B. 
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FIG. 1. Complex activity plane for a lattice or contin- 
uum fluid system at temperatures T > Tc, above gas- liquid 
criticality. The Yang-Lee edge singularities are located at 
z = Air) ± iz'AT) while the repulsive-core singularity lies 
on the real negative activity axis z = zq{T). (Note that 
the branch cuts running from the Yang-Lee edges will not, 
in general, be linear as shown, purely schematically, here: in- 
deed, for a simple lattice-gas the cuts lie on circles centered 
at z ^ 0.) 



FIG. 2. Plot of t/je""'' (bold curve) versus tp to elucidate 
solutions of the saddle-point equation (18) in the case of only 
single-site repulsive hardcores with attractive potentials 
of strength vo (= —Uo) for various temperatures, l3c and 
132 and positive activities z. The open circle marks the point 
of inflection which serves to locate the (classical) gas-liquid 
critical point; the dashed line is the associated tangent. The 
solution lines have slopes —z in accord with (18). 
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FIG. 3. Plot of xlie~^ to illuminate the saddle-point equa- 
tion (18) as in Fig. 2; but note the change in scales. The dom- 
inant repulsive case with i7o = wo — 'So > is illustrated for 
positive activities z [lines (a) and (b)] and negative activities: 
(c) and (d) and the tangent (dashed line) which locates the 
repulsive core-singularity at z — zo{T) < with field ^o{T). 
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